Spectral Decomposition of the Tent Map 
with Varying Height 

Suresh Subbiah 1 and Dean J. Driebe 1 ' 2 
1 Center for Studies in Statistical Mechanics and Complex Systems 
The University of Texas at Austin 
Austin, Texas 78712 USA 
2 International Solvay Institutes for Physics and Chemistry 
1050 Brussels, Belgium 

Abstract 

The generalized spectral decomposition of the Frobenius-Perron operator of the 
tent map with varying height is determined at the band-splitting points. The de- 
composition includes both decay onto the attracting set and the approach to the 
asymptotically periodic state on the attractor. Explicit compact expressions for the 
polynomial eigenstates are obtained using algebraic techniques. 

1 Introduction 

A natural description of the dynamics of chaotic systems is in terms of evolving probability 
densities On this level the time evolution in maps is governed by the linear Frobenius- 
Perron operator and the dynamical problem is solved by the determination of the spectral 
decomposition of this operator. In recent years, several authors have constructed complete 
and explicit spectral decompositions of the Frobenius-Perron operator of a variety of model 
systems [||, |j. The most useful decompositions contain in their spectrum the decay rates 
characterizing the approach to equilibrium of the system. For one-dimensional piecewise- 
linear Markov maps such decompositions are constructed in function spaces spanned by 
polynomials. The dual space of these polynomials are generalized function spaces and so 
the decompositions are known as generalized spectral decompositions 
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Chaotic systems often contain a control parameter that characterizes the strength of 
the chaos. A simple model system with such a control parameter is the well-known tent 
map with varying height p|. The tent map with height h on the unit interval is given by 



T(x) 




< x < 



I) 



- x) I < x < 1, 

where the parameter a = 2h. In Figure 1 the map (|l]) with height y/2/2 is shown. Note 
that the map acts on the unit interval [0,1) but has images on [0,h]. As characterized 




II 



III 



IV 



Figure 1: The tent map at a = For this value of a iterates of the critical 
point define four intervals around which the dynamics is organized, as discussed in 
Section 2. 



by its Lyapunov exponent, log a, the map ([!]) switches abruptly chaotic as the height is 
raised past 1/2. We are interested in the chaotic regime where 1/2 < h < 1, i.e., 1 < 
a < 2. Similar to the well-known universality of the quadratic map, it has recently been 
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reported || that the tent map also governs the low-dimensional behavior of a wide class 
of nonlinear phenomena. Specifically, it was found that the dynamics of the Ginzburg- 
Landau equation, in its description of the modulational instability of a wave train, is 
reducible to the tent map. Under this reduction, varying the height in (|l|) corresponds to 
varying the wavelength of the initial modulational instability. 

Our interest is in the statistical properties of the iterates of the tent map and in evolving 
probability densities. The Frobenius-Perron operator, U, corresponding to a map, S(x), 
defined on the unit interval evolves a probability density, p(x, t) by one time step as 



p(x,t + 1) = U p(x,t) = / dx 5(x — S(x')) p(x' , t). 

Jo 



(2) 



Evaluating the integral gives a sum of contributions from the inverse branches of S(x). 
The Frobenius-Perron operator corresponding to the map (jl|) acts explicitly on a density 
p(x) as 

'~e(%-x), (3) 



Utp(x) = - 
a 



x\ i a — x 
P\ ~ )+P[ 



where 



The step function appears here because the map has no inverse images for x > a/2. For 
a = 2 the map has images on the whole unit interval. For this value of a the invariant 
density (being the stationary solution of (|3|)) is uniform on the whole unit interval. As 
a is lowered the invariant density is supported only on a subset of the interval [a(l — 
a/2), a/2], i.e., from the second iterate to the first iterate of the critical point x c = 1/2. 
The invariant density is discontinuous at all values of the trajectory of the critical point. 
If the critical trajectory is periodic (or eventually periodic) there will be a finite number 
of discontinuities. 

For a > \/2 the invariant density has nonvanishing support on all of [a(l — a/2), a/2]. 
As a is decreased past this critical value, a%, the invariant density breaks up into two bands 
with a gap in the middle. Decreasing a past a<i = 2 1 / 4 causes the invariant density to 
break up into 4 bands. In general the invariant density has 2 n bands as a is decreased past 
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These values of a are called the band-splitting points M. This is illustrated 



in Figure 2. At the band-splitting points the critical trajectory is eventually periodic 








(b) 


T (XcK 












//t (5) (x c) \ , 








//^-~~~ (3T-\ \ //, 








„i 4 > / / / 




T (Xji/ / / 
— T (Xc>/ / 




(8) / 








\T (x c ) 



Figure 2: Bifurcation plot of the tent map with varying height showing the for- 
mation of bands as the height, parametrized by a = 2h on the horizontal axis, is 
lowered. The vertical axis is the subset [a(l — a/2), a/2] of the one-dimensional 
phase space. In (a) the shaded regions indicate where the invariant density has 
nonvanishing support. At this resolution we see up to the formation of the four 
bands at a = 2 1 / 4 but higher band-splitting points are not resolved. In (b) iterates 
of the critical trajectory as a function of a are plotted. These are seen to determine 
the band structure in (a). 



and the invariant density is constant in each band and thus piecewise-constant over the 
unit interval. But a general initial density will also have persistent oscillating components 
among the bands. This feature is known as asymptotic periodicity Q. 

We want to determine the spectral decomposition of the Frobenius-Perron operator 
so that we may expand a density or correlation function in terms of its eigenmodes. In 
order to do this we need the dual states or left eigenstates of U. These correspond to the 
right eigenstates of the adjoint of U, which is known as the Koopman operator The 
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Koopman operator, K = W, acts on a phase space function A(x) as 

KA{x) = A(T(x)), (5) 

where T(x) is the rule for the map, such as (JlJ) . For decompositions of one-dimensional, 
chaotic, Markov maps in spaces spanned by polynomials the Koopman operator has eigen- 
states that are generalized functions or eigenfunctionals || ||, f|] . 

The decay of the x-autocorrelation function for the tent map at the band-splitting 
points and at values of a close to band-splitting points was calculated in ||. The tent map 
has also been studied by Dorfle Q at arbitrary values of a in several function spaces but he 
does not provide explicit complete spectral decompositions. The asymptotic periodicity of 
the system has been studied by Provatas and Mackey [fTCfl . Since as t — > oo the density is 
supported only within [a(l — a/2), a/2] all these authors have only considered the map in 
that region and neglected transient behavior onto it. That is sufficient if one only considers 
the behavior of time correlation functions; but for the general evolution of densities and 
observables in a nonequilibrium statistical mechanics context this transient behavior must 
not be neglected since, as will be seen, the slowest decay modes originate from this part 
of the dynamics. 

In the next section we construct the spectral decomposition of the Frobenius-Perron 
operator at the first band-splitting point. In Section 3 the decomposition at all the band- 
splitting points is constructed using the self-similarity of the map at higher band-splitting 
points to the map at lower band-splitting points. 

2 The first band-splitting point 

From (|l|) the map at the first band-splitting point corresponding to the height h = \[2/2 
or a = \[2 is 

f V2x < x < 4 

Ti(aO = <^ ^ 2 (6) 

1 V2(l-x) \ <x< 1, 

which is shown in Figure 1. The dynamics is organized around four intervals determined by 

the trajectory of the critical point. At the third iteration the critical trajectory settles onto 
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the fixed point, x* = 2 - y/2. The four intervals: I = [0, tJ 2) (x c )), II = [if \x c ), Tj 3) (x c )) 
III = [Tj (x c ), (x c )) and IV = [T^(x c ), 1), define a minimal Markov partition for the 
map. (These intervals are indicated in Figure 1.) Any point in the interior of interval IV 
is mapped onto some point in interval I in one iteration. Under successive iterations all 
points in the interior of I are eventually mapped into II. Any point in interval II maps 
into interval III in one iteration and any point in interval III maps into II in one iteration. 
Thus the union of the intervals II and III form the attracting set f2. 
The Frobenius-Perron operator for the map (jpT) acts on a density as 



x \ / v I — x 



p\-7n)+p\ 



fV2- 



/V2 



G(^--x). (7) 



V2J r V V2 

A general initial density continuous over the unit interval develops discontinuities at the 
endpoints of the four intervals described above. We thus choose a function space to 
consider Ut 1 in as the space of piecewise-polynomial functions where the pieces are the 
four intervals described above. The invariant density has support only in $7. 

The fact that points oscillate between the intervals II and III means that a general 
density will have a persistent oscillating component between these two intervals under 
time evolution. This property we will sometimes refer to in the present context as the 
"flip property" , since the part of the density in II will all be in III (with stretching) in the 
next time step and vice-versa. The parts of the initial density with support in the intervals 
I and IV will decay onto f2. Since we will use the eigenfunctions on f2 to determine those 
on its complement, we first consider the decomposition on 0. 

2.1 Evolution on the attractor 

For convenience f2 is stretched onto the interval [0, 1). At the end of the computation we 
will rescale all results back to £1. The linear function that makes the stretch is 

<p(x) = (2/x*)x-V2, (8) 

where T^\x c ) < x < T^(x c ). The transformation @ is a homeomorphism so that it 
begets a new map R4 topologically conjugate to the part of Tj on O as R4 = 4> Ti o (ft^ 1 
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Ri (x) 



(9) 



given by 

j V2x + x* 0<x<x*/2 
\ y/2(l-x) x*/2 < x < 1. 

where x* = 2 — \[2 is the fixed point of the map Ri(x), which is the same as the fixed point 
of the map Ti(x). Under this transformation, the intervals II and III are stretched to the 
intervals A = [0,x*) and B = [x*, 1) respectively. The map Ri(x) is shown in Figure 3. 

1 



Rl(x) 




Figure 3: The rescaled Tent map at the first band-splitting point. 



The Frobenius-Perron operator corresponding to the rescaled map Ri acts on a density 



as 



Ur^x) 



1 

71 



p(^) +P ( x -^f)e(x-x* 



V2 



(10) 



The flip property of Ti is inherited by Ri in that the inverse image of A is B and vice- 
versa. This suggests that a simpler analysis will be obtained by considering the map 
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corresponding to two iterations of Ri. This map, Gi = Ri o R 1; is given by 

i-2x + x* < x < x*/2 

2x-x* x*/2 < x < (1 + x*)/2 (11) 

-2x + (2 + x*) (l + x*)/2 <x < 1. 

The flip property of Ri means that Gi is metrically decomposable into two independent 
maps on the intervals A and B, as is clear from Figure 4. From now on in this section 
we shall drop the subscript 1 on the maps Ri and Gi; it being understood that we are 
referring to these maps at the first band-splitting point. 




Figure 4: The map Gi = Ri o is metrically decomposable into two parts each 
conjugate to the tent map with unit height. 

The map G restricted to A is just a rescaling (with a flip) of the tent map with full 
height, To, i.e., the map (|]) with a = 2. This is expressed in terms of a topological 
conjugacy as 

G A (x) = (/) A 1 (x)oT o (x)o0 A (x), (12) 
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where the conjugating function 4>\(x) is 

cj> A ( x ) = l-^-, ( 13 ) 
x* 

and x E [0, x*). Similarly the map on B is topologically conjugate to To as 

Gb(x) = ^ b 1 (x)oT (x)o^( 2; ), (14) 

where the conjugating function 0b (x) is 

<fa(x) = (V2/x*)x- y/2, (15) 

and here x £ [x*,l). These conjugacies are useful for us because the spectral decom- 
positions of maps that are topologically conjugate are simply related, as is reviewed in 
Appendix A. The generalized spectral decomposition of Tq has been previously deter- 



mined [11, 12] and is reviewed in Appendix B. 

Following the discussion in those appendices gives the right eigenvectors of Ga and 
Gb as 

r*)o A = > 2 ,(^)*A (16a) 
I<Wg a = ^E 2j+1 (^) X A, (16b) 



I 2 ^')gb = : §B 2J (^ 7 X) XB (17a) 



x* 3 V V2x* 

|02i+l>G B = ^ E ^+ l {^^ X ~ X *^j^ ^ 17b ) 

where the associated eigenvalue is the argument of the ket vector with |02j+i) meaning 
a null eigenpolynomial of degree 2j + 1 and xa an d Xb are indicator functions on the 
intervals A and B respectively. Due to the metric decomposability of G these states are 
eigenstates of Uq as well. 
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Similarly, the left eigenvectors of Ga and Gb are the generalized functions 

<2-*|G A - ^^ [^(x-O-ffl^W] (18a) 
(0 2j+ i|g a = (18b) 



(2- 2j |g b = ^ ^-[^(x-lJ-^Cx-x*)] (19a) 

"I f±*\* +2 M + i) 



(0 2j+1 |c B = (1%) 

where the definitions of <5± are given in Appendix B. 

Since C/q = £7r the spectrum of C/r is a subset of {0, ±2 -J '}. Consider a non-zero 
eigenvalue, 2~ 2j , of C/q- There are two eigenvectors associated with this eigenvalue, each a 
polynomial of order 2j. Since the function space on which C/r acts has two basis elements 
for each degree j, i.e., an j th degree polynomial in A and an j th degree polynomial in B, 
there should be either two eigenvectors or one eigenvector and one Jordan vector that are 
polynomials of degree 2j, associated with either one or both the eigenvalues {+2~ ] , — 2 _J }. 
Since Uq does not have any Jordan vectors it follows that C/r doesn't either (for non-zero 
eigenvalues). The eigenvalues of C/r cannot be twofold degenerate since that would imply 
that all the eigenvectors of Uq are also eigenvectors of C/r, which is impossible since G is 
metrically decomposable and R has the flip property. Therefore the non-zero eigenvalues 
of R are +2~' and 

The eigenvectors of C/r with eigenvalue ±2~- J are in the eigenspace spanned by the two 
eigenvectors of Uq corresponding to +2~ 2j \ Thus they will be linear combinations as 

l+2- J >R = \ (|+2- 2j )g a + c 3 \+2^) Gb ) (20a) 

and 

|-2-*>R = \ (|+2- 2 ^) Ga +dj \+2-^) GB ) , (20b) 
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where the coefficent of 1/2 is put for convenient normalization. To determine Cj we use 
that |+2 — is an eigenvector of Ur as 



U R |+2-^') Ga + Cj |+2- 



-2j\ 



2~ 3 



+2~ 2j ) GA + Cj \+2 



-2j\ 



The flip property tells us that 



Ur 



+2 



-2r 



g a - 



(21) 



(22) 



Substituting the explicit form ( |16a|) of |+2 2j )g a an d (|17a|) of |+2 2:, )gb i R (H) an d 



solving for c,- we find that Cj = 2 J . A similar analysis shows that dj = —2 3 . Hence 



l±2-') R = JL ft, 



v 7 ^ (x-x* 



XB 



The invariant state, corresponding to the invariant density of Ur is 

1 



I+i)r 



2x- 



■(xa + V2~Xb) 



(23) 



(24) 



This state carries all the probability under evolution of Ur and any density will have this 
component. The state 



-i: 



R 



2x' 



-(XA-V2 



XB) 



(25) 



is the asymptotically periodic state. Only it and the invariant density survives as t — > oo; 
but | — 1)r, like decaying states, doesn't carry any probability. It does keep the memory 
though of the projection of the initial density on | — 1)r, which is a special property of 
asymptotically periodic systems ||] . 

Now we consider the null space of Ur. The map G has two independent null vectors 
(one in A and one in B) for each odd degree. This implies that R can have either a 
corresponding 2x2 Jordan block or have 2 independent eigenvectors for each odd degree. 
The latter case is not possible since null vectors of R cannot have support in interval B 
because only one of the terms on the rhs of (|io| ) acts on functions in B. Thus there is a 
2x2 Jordan block for each odd degree associated with eigenvalue 0. 



11 



Consider the action of Uq = U R on a null state, |02j+i)ga> °f Ga as 



Ur U r \0 2j+1 ) Ga 



0. 



(26) 



The function inside the square brackets has support only in B, and Ur acting on any non- 
zero function with support in B cannot vanish in one iteration. Thus |02j+i)r = |02j+i)ga 
is a null vector of Ur with explicit form given in (|16b| ) 

The Jordan vector, |0j 23+1 )r, associated with this eigenvector satisfies 



Ur\0j 2j + 1 )r = |0 2 j + l)R. 

We may choose the Jordan vector to have support only in B as 



(27) 



(28) 



where f]2j+i is a constant to be determined. To determine ??2j+i we apply Ur on (|28|) and 
use (27) and the explicit forms ( |16b| ) and ( |17b| ) (remembering that |02j+i)R = |02j+i)G A ) 
to obtain rj2j+i = — 1. 

2.1.1 Left eigenstates of Ur 

The left eigenstates of Ur with non-zero eigenvalues may be determined by an approach 
similar to the one used to find the right eigenstates. The result is 

(2x*) 2j 



(±2- 



\r 



m 



fpy-i) 



(x) -5 { y~ 1] (x-x*) 



±5 



(2j-l) 



(x-x*)^^" 1 ^^-! 



(29) 



The left states with zero eigenvalues are like the right states identical to those of Uq 
(except for the factor of — 1 for the dual state of the Jordan vector) as 



(°2?+i|r 

(0j 2j + 1 |R 



(°2j+llG A 
~ (°2j+iIg b • 



(30a) 
(30b) 



Note that (0j 2j+1 |r is an eigenstate of the Koopman operator and (02j+i| R is a Jordan 
state. 
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2.2 Decay onto the attractor 



As noted before, initial densities with support in the intervals I and/or IV will decay into 
17. Consider a density with support only in IV at t = 0. At t = 1 the density has support 
only in I. Thus any eigenvector of Ut 1 with support in IV can only have eigenvalue 0. To 
determine such an eigenvector we write an ansatz for it as 



\0j)w = fi,j(x)xi + fu,j( x )xn + fm,j{x)xin + fw,j(%)xiv, 



(31) 



where the subscript IV on the ket denotes that it describes decay out of interval IV and 
we take fwj(x) as a polynomial of order j. Applying Ut^ to (|3l|) and collecting terms 
that are multiplied by the same indicator function we get 



U Tl \0j)w 



+ 
+ 



x 



x 

X 



XII 
XIII- 



(32) 



Since |0,-)rv is a null eigenstate the coefficients of each of the indicator functions must 
vanish. Since /rv j is a polynomial of degree j, it follows that f\t and fnij are also 
polynomials of degree j. We choose fnj to be zero. Clearly, fmj = /rv,j and choosing it 
as (x — l) 3 fixes fjj. We thus have 







'i/iv 



-l) j+1 x j X i + (x- l) j (xm + Xiv)- 



(33) 



We note that because of the degeneracy associated with eigenvalue in interval IV, this 
choice of eigenvectors is not unique. The left eigenstates (given below) associated with this 
choice of right eigenvectors is unique and therefore when we expand an arbitrary density 
in terms of the right eigenvectors, the expansion coefficents are uniquely defined. 
The action of Ut 1 on a function with support only on I is given by 



UtMWxi] = -^f(x/V2)( X i + xn) 



(34) 
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Acting with Ut 1 on a monomial in I gives 

U T1 WX1] = T^-^(XI + Xli), (35) 

so that the eigenvectors of Ut 1 with support in I are monomials in I. To determine their 
form in the other intervals we again write an ansatz as we did for the eigenvectors with 
support in IV as 

|2- (j+1)/2 )i = x> X l + 9uA x )xn + 9iiiA x )xm, (36) 



where the associated eigenvalue, appearing as the argument of the ket, is seen from (35). 
In equation ( |36[ ) and below the subscript I on a ket implies that it describes decay out 
of region I. Note that the j = mode here is the slowest decay mode in this system. 
Since Ut 1 does not raise the degree of the polynomial it acts on, giij(x) and gmj(x) 
are polynomials of degree j. Applying Ut 1 to ( j36|) and using ( |35|) and the fact that the 
function on the rhs of ( |36| ) is an eigenvector with eigenvalue 2~ v+ 1 )/ 2 ) we find that 

UT 1 (gu,j(x)xn + gm,j{x)xm) = 2~ (j+1)/2 (g U j(x)xn + gm,j{x)xm - x j Xli)- (37) 

The determination of guj(x) and ginj(x) is described in Appendix C 
The explicit form of the first few of these eigenvectors is 

|2 " /2)i = (38a) 

I2- 1 ): = xxi- ^11)0 + ^|-l>n + ^|0>n (38b) 
M h - x Xi j^— |l)nH ^ 1 — 1/n H ^ |U)n (38cj 

2 — 12 )n + ^iT7I) 1 h - (38d) 

The subscript f2 on a ket implies that it is a right state of C/r which has been rescaled into 
the attractor, f2, of Ut 1 . The explicit form of these states is given in Table 1. 
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2.3 Left eigenstates of 



It is easily seen that the left states given by 



1) 



(39) 



form a bi-orthonormal set with the null states in ([33]) . These left states are also orthogonal 
to all the other right states of Ut 1 . The left states associated with the transient decay 
states out of I are given by 



<2- {j+1)/2 |i = \\( 



-l) j 6 i i ) (x) + 5 U) i 



x-1) 



(40) 



The left states of J7r are not the left states of Ut 1 , even though the right states of ?7r 
are also right states of Ut 1 . This is because Ut 1 acting on a density contained within 0, 
will continue to be in f2. But in general the Koopman operator, K^, acting on a function 
with support only in 0, will result in the function having support outside O too. The left 
states of Ur scaled back to form a bi-orthonormal set with the right states contained in 
f2, but they are not orthogonal to the transient states that decay into f2. To make them 
so we use Gram-Schmidt orthogonalization. The results are given in Table 1. 

2.4 The spectral decomposition 

Using all the eigenstates and eigenvalues given in Table 1 we may write the action of 
in terms of its spectral decomposition as 



oo 

|i)<i| n + (-i)*|-i)(-i| n + E 0r<* +i v 2 ) \ 2 -mm ){2 



-(2j+l)/2 



3=0 



3=1 



| + 2-J)(+2-^| n + |2-^)(2-^|i 



+ E(-2- J ) t |-2- i )(-2^b + E^I°2 J +i)(0j 2j+1 b 

j=l j=0 

oo 

+ E^* [|0 2 i+i)(0 2j +ib + |0j 2 , +1 )(0j 2j+1 b + \0j)(0j\w 

3=0 



(41) 



where the subscript on the bra states also identifies their dual ket states. 
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eigenvalue 



degeneracy 



symbol 



eigenvector 



1 



1 



H>n 
(l|n 



XII 



+ V2 X 



in 



(xii + xm) 



E2°=o (a , fc (2-( fc+1 )/ 2 | I + a' 0ifc (0 fc | IV 



|-l>n 
<-l|n 



xii - V%x 



in 



(xii - Xlll) 



Er=o (b M2- {k+1)/2 \i-a'^ k (0 k \ w 



|+2- 
(+2 



—J 



(^) 4j - 2 r 



?(2i-i) 



(i-TWy-^^i-i*) 



Er=2,(^(2- (fc+1)/2 |i + a;- fe (0 fc 



IV 



2-> 



<-2- J '|n 



^ (^) xii - #£ 2 , (^) XIII 



(2j)! 



?(2i-l) 



(x-T( 2 )(x c ))-5l 2i-1) (x-x*) 



15 (x - x*) + 5i 2j - 1} (x - T( 1 )(x c )) 

Er= 2 ,(^(2- (fc+1)/2 ii-4 fc (o fc | IV 



OC 



(2x2 
Jordan 
blocks) 



l°2j+l)n 
(°2i+iln 



Xii 



|0j : 



aj+i / 



-Er= 2 , + i^(2- (fc+i)/2 ii 



ii 

x* 



XIII 



0./: 



2j + l 



(**) 4j+2 1 A^+ 1 )( T tW^t ^ 



E 



fc=2j+l Cj, 



| 2 -(J+ 1 )/2) I 
( 2 -0'+i)/2| I 



? Xi + E i=0 a i,P~ % ) + M " 2_i ) + c iM) 



i ((_i)i^(i)( x ) i) 



|0i)iv 
(OillV 



(-iy+v X i + (x-iy(xiii + xiv) 



Table 1: Elements of the spectral decomposition of the tent map at the first band splitting point. 
The constants a itj = (2- i |n(2)-0'+ 1 )/ 2 ) Il b itj = (-2-'| n (2);«+ 1 )/ 2 ) I and * d = (0i|n(2)-W +1 >/ a )i 
and are given in Appendix C. The constants = (+2~*|n0j)rv, = (— 2 ~ l |n0j)rv = ~ a 'i,j 
and j' t j = (Oj; |n0j)iv and <p(x) is defined in @. 
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3 Higher band-splitting points 



We can determine the spectral decomposition of the tent map at any band-splitting point 
(bsp) by generalizing the approach used in the previous section for the first bsp. For 
the rescaled map at a = V2 we found that by considering its square it separated into 
two parts that were directly related by a simple change of scale (including a reflection 
for one part) to the tent map at full height where a = 2. This relationship also holds 
generally [0| between the map at the n th bsp, (a n = 2 2 ), and the (n — l) th bsp, 
(a n _! = 2 2 ~ ( - 1) = a n ). 

The tent map at a n on the interval [Tn\x c ),Tn\x c )), which contains the attractor, is 
first stretched to the unit interval [0, 1) to make the rescaled map R n . The linear function 
that makes this stretch is 

i / \ 2x — a n (2 — a n ) . . 

<Pn{x) = 7 tt — • (42) 

«n(«n - 1) 

Under ([42^) the map T n transforms to R n = (fr n o T n o and is given by 

{2 — a n (l — x) < x < (a n — l)/a n 
K '- 11 (43) 

a n (l-x) (a n - l)/a n < x < 1. 

We then compose R n with itself to obtain G n = R n o R„. As is illustrated in Figure 5 
for the 2 nd bsp, It can be shown easily that in general: 

(a) G n (x) in the interval X n ^ = [0,Rn 4 ^(x c )) is topologically conjugate to R n -i (the 
rescaled map at higher height) in the interval [0, 1) as 

Rn-l = 0n,A o G njA ° 4>n\i ( 44 ) 

where 

x 

4>n.A = 1 1 rr, X G X n .A- (45) 

a n {a n - I) 

(3) 

(b) G n (x) in the interval X n ^ = [Rn (x c ),l) is topologically conjugate to R n _i in the 
interval [0, 1) as 

Rn-l = 4>n,B o G„ iB ° <j>nB: ( 46 ) 
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where 



x - (2 - aj 

VB = ^ ) x & A n,B- 



a,, 



(47) 



G 2 (x) 




Figure 5: At the second band-splitting point, G2 is conjugate in X2,a and X23 to 
Ri shown in Figure 3. The central transient region is shown as A^c- 



At the n th bsp the critical trajectory is eventually periodic with period 2 n ~ l (n > 1), 
and therefore the number of discontinuities that appear under time evolution of an initially 

imply that the band structure at the n th 



smooth density is finite. Equations (44) - 
bsp consists of 2 scaled copies of the band structure at the (n — l) th bsp separated by 



the interval [Kn\x c ), Rn 0j (x c )) = X n> c, which makes one band. Therefore we have the 



,(3) 
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following recursion relation for S n , the number of bands at the n bsp for the map R ra , 

S n = 2S n -! + l, (48) 

for n > 2, where S\ = 2. Solving this recursion relation gives 

S n = 2 n + 2 n - 1 - 1, (49) 

for n > 1. Of the S n bands, the invariant density has support only on 2 n bands. The 
density is transient on the other 2 n ~ 1 — 1 bands. The function space we consider at the 
n th bsp is piecewise polynomial, with each piece extending over one band. 

3.1 Decomposition on the attractor 

Associated with the 2 n bands that make up Q n , the attractor at the n th bsp, are 2 n 
eigen/ Jordan vectors of each polynomial degree j. We denote an element in the spectrum 
of Uji n associated with states on Q n as X% j and the right eigen/Jordan vector associated 
with it as |AJJ ■). The superscript n stands for the order of the band-splitting point. 
The index k is an integer from 1 to 2 n and distinguishes between the 2 n independent 
eigen/Jordan vectors of degree j, at the n th bsp. 

Since G„ 5 a and G„ 5 b are topologically conjugate to R n _i all three share the same 
spectrum. The map G ra is the union of G nj A ; G Ki b and the part of G n on the central 
transient interval X n £. The spectrum of R ra on £l n is determined from the spectrum of 
G n on £l n . Consider a non-zero eigenvalue of R n -i (we will show later that as at 

the first bsp the right states corresponding to non-zero eigenvalues on the attractor are 
even-order polynomials), which is also an eigenvalue of G n with degeneracy 2 (one for 
G n; A and one for G n> B)- Using arguments parallel to those that we used at the first bsp 
we deduce that R n has two distinct eigenvalues AjJ^j and — ^/a^J. By induction 
R n+ i on Q n+1 has in its spectrum the eigenvalues -(A^j) 1 / 4 , + i (.K72j) 1/4 

and — ^(A^j) 1 ^ 4 - Thus from the non-zero eigenvalues at the th bsp we can determine the 
non-zero eigenvalues at the n th bsp by taking square roots of 2 -2 - 7 recursively, n times. 
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Thus 

*bi = (2 J exp (-*-) ' (50) 

where k = 1, 2, . . . , 2 n and j = 0, 1, 2, . . . . Note that for k < 2"- 1 

Afc,2j = ~^k+2™~ 1 ,2j- (51) 

At the th bsp the eigenvalue is infinitely degenerate with an associated eigenpoly- 
nomial of odd degree for each occurence of the eigenvalue. In the previous section we saw 
that at the 1 st bsp there are an infinite number of 2 x 2 Jordan blocks associated with 
eigenvalue zero, one each for each odd degree. We will show by explicit construction below 
that this trend continues and at the n th bsp, R n has a 2" x 2 n Jordan block for every odd 
degree. 

3.1.1 Right states on the attractor 

To determine the right states on Q n with non-zero eigenvalues we begin similar to (20) by 
writing the right states at the n th bsp in terms of those at the (n — l) th bsp as 

\K,2j) = l A fc,2i)G n , A + c fe,2jl A fc,2i)G n , B ( 52a ) 
l^fc+2™-\2j) = l^fe,2j)G n , A + d k,2j\K,2j)G n , B , (52b) 

where k = 1, 2, ... , 2 n_1 . In equation (52) and in the remaining part of this subsection 
kets without any subscript denote right states of R at the appropriate bsp. Conjugacies 
( fn| ) and ( |46| ) imply that the eigenvectors of G nj A and G n) B are related to those of R n _i 



as 



|A£,2i>G„, A = U rA \\ n k ^) (53a) 



Using this in (52) gives 



|Afo) = U^\\Z))+cl 2 p^\\-)) (54a) 
\K + 2^,2j) = U^X^ + dfaU^Xffi). (54b) 
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The state U 6 -i |AJ? 2 j) nas support only in X n ^ and U^-i | A£ 2 j) has support only in X Ui b- 



We now use the relation 



BW^B = ^.A> ( 55 ) 



which implies that 



^JVl^V (56) 
(Even though (|55|) and (|56| ) are written specifically at the band-splitting points, they are 
valid for all a £ (1,2].) Acting on (54) by Un n and using (^) we get 

AfolAfo) = U^U^X^+dfaUf-iJXfe) (57a) 

where (|5l| ) has been used on the lhs of ( |57b| ). Since Un n has the flip property that any 
function with support only in Al H) b will go entirely to X n ^ in one iteration and vice- versa 
we know that Un n U,-i jA^l) has support only in X Uy B and U,-i |AJJ g 1 ) has support only 
in X„ a- Multiplying ( |54a| ) by X^ 2 j and fl54h|) by — A^ 2 j an d identifying the components 
with support only in X n) x with the corresponding components in (57) gives 



c k,2j U <j,-^\ X k,2j) = X k,23 U ^-\\ X k,2j) ( 58a ) 

dh^WU) = -A^.^JA-j), (58b) 

showing that c% 2 j = A)? 2 j an d d k2j = — A^ 2 j- Using this result gives the pair of recursion 
relations 

Mm) = K-i + ^Va) 1 ^^ (59a) 

|A^n- 1)2i ) = (U^ - AjfoE^) |A-j), (59b) 

which express the right states at the n th bsp in terms of those at the (n — l) th bsp. 

These recursion relations can be solved to write the right states at the n th bsp in terms 
of the right states at the th (a = 2) bsp. For notational convenience we define 

B, = U+£ (60) 
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of Aj's and B^ 2j - j's, where if Oi = 1 then the i th location in the n-product (counting from 
the right) will be taken by B£ 2i and if <7j = then the i th location will be taken by A*. 



Let G{ (i = 1, 2, . . . , n) be either or 1 and we define n o - n0 - n _ 1 ... - 1 to be an ordered n-product 

\n J 
5 fc,2j,i 

r f"\ Irfwi Kir 

Solving (59) gives 

\^k,2j) = y^IIgnO-n-l-o-l 1^1.27) ■ 
M 

The sum in ( |6T|) is over all possible cr-strings of O's and l's of length n and so consists 
of 2 n terms (n-products). The order of the operators in each n-product must be strictly 
observed since the operators involved do not commute. 

To illustrate (^) we write it out explicitly for n = 1 and 2. For n = 1 ( |6l| ) gives 

\^k,2j) = Al|A° 2< j) + ^l,2j,l\^l,2j)- (62) 

This agrees with the expression (^3|) (corresponding to k = 1 and = 2) we had for the 
right eigenstates at the first bsp. For n = 2 (^3l]) gives 

l^fc,2j) = A 2 Ai|A5 2j -) + A 2 B| ]2 j,i|A ( i , i 2j) 

+ B l,2i,2 A l|Al,2j) + B fc,2j,2 B fc,2j,l|Al,2j)- (63) 

Now we prove by induction that there is a 2" x 2 n Jordan block associated with the 
eigenvalue at the n th bsp for each odd order 2j + 1. In section 2 it was shown that this 
statement is true for the first bsp. Assume that this statement is true at the (n — l) th 
bsp. We denote the Jordan vectors as jO^j+i) where k = 2,3, ... ,2 n_1 (IO^j+i) is the 
eigenvector of the block). They satisfy 

UrM^ +1 ) = \O n k Zl >2j+1 ), k + \ (64a) 
^RnlO^i+i) = 0. (64b) 

Since we are assuming that C/r„_ 1 has 2 n_1 x 2 n_1 Jordan blocks, the conjugacies (p4|) and 
( |4"6| ) imply that G ni A and G nj B both have 2 n x 2 n_1 Jordan blocks with states given by 

|0fc, 2 j+i>G n , A = AnlO^+i) (65a) 
K,2 j+ l)G n , B = B n |O0j +1 ), (65b) 
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where k = 1, 2, ... , 2 n_1 . Since G n on Q n has two 2 n_1 x 2 n_1 Jordan blocks for each j, 
R n on fi n can either have two 2 n_1 x 2 n_1 Jordan blocks for each j or have one 2 n x 2 n 
Jordan block for each j. The first case implies that J7r„ have two null vectors |0™ 2 j+i)G n A 
and |0™ 2 j+i)G n B f° r eacn i- But £^R n 1 2j'4-l ) G n b ^ since no function with support in 
X Ut B can vanish in one iteration under U^ n . Therefore Un n has a 2" x 2" Jordan block 
for each odd degree 2j + 1. This completes the proof by induction. 

A null state of C/r^ has to be a null state of Uq u too. Therefore the null state of U^ n 
for each j is given by 



We use the relation 



which implies that 



0?, 2i+ i>= A -l°i,2i + l>- (66) 



R « </VA = 0n,B R n-1 ( 67 ) 



U Rn A n = B n U Rn _ 1 . (68) 

Unlike d56[), equation ( |68|) is true only at the band-splitting points. Equation (^) can be 
used to verify that if we act on both sides of (|66|) by Un n its rhs reduces to zero. 

The right state |0™ 2 j+i)G n B ^ s a good candidate for the Jordan state lO^j+i) since 
under one iteration by C/R n it will have support only in X n> A and under two iterations of 
Un it will vanish. Tentatively, from (|65b|) we write 



|0^ 2 ,+i) = B n |0^ 1 +1 ). (69) 

This guess for the Jordan state can be verified by using relation ([5(]). Similarly ( p8| ) can 
be used to show that 

|0Wi)= 1 n|0^ 1 +1 ). (70) 



and equation (56) can be used to show that 

-*4,2j+l) = LJ n\v 2t 2j+l/- 



0Wl)= B n|02.2, 1 +l)- (71) 
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In general, we find 

K,2 j+ i) = K\0^/\l2 j+ l) for A; odd (72) 
K2 j+ i) = B„|0^ 2i+1 ) for k even 

where k = 1, 2, . . . , 2 n and [g] is the ceiling function (q if it is an integer or else the next 
greatest integer). These recursions are illustrated in Figure 6 up to n = 2. This recursion 
relation can then be solved to write the eigen/Jordan vectors at the n th bsp in terms of the 
null vectors at the th bsp. To write down a compact solution we define n o - ji(Tn _ li .. CT1 , which 
is similar to the n CT „o- n _i...<T 1 previously defined. We define n (Tn0 - n _ 1 ... (Jl to be an ordered 
n-product of Aj's and Bj's. If = 1 then the i th location in the n-product (counting from 
the right) will be taken by Bj and if <7j = then the i th location will be taken by Aj. With 
each H an <7 n _ 1 ...<j 1 we associate a binary number formed from the string of l's and 0's as 
k = a n a n -i . . . a\ + 1. Solving the recursion relation (|72|) we get 

|0re,2j+l) = no-„cr n „ 1 ...cri|0 ( j ) j 2j+l)) (73) 

where |05 2j+i) ls the null vector of degree 2j + 1 of the tent map with full height, and k 
here is the decimal equivalent of the binary k, which ranges from 1 to 2 n . 
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3.1.2 Left states on the attractor 



We obtain the left states, (A£ -|, which are orthonormal to the right states given by ( pl| ) 
and ( ff3D by taking the duals of those expressions. The dual expression of (^Tj) gives the 
left states corresponding to the non-zero eigenvalues as 




where the n-product here is of [A i )t's and ((BjJ 2ji ) 1 )^'s and the factor of l/2 n is put 



for normalization. The left states corresponding to the Jordan vectors associated with the 
eigenvalue are 



3.2 Decay onto the attractor of the rescaled map 

We saw in section 2.1 that Ri has no transient bands. The band structure of G2 consists of 
two scaled copies of that of Ri separated by the central interval -X^c- This central interval 
is a transient band of R2. Since G3 5 a and G3 ; b have a band structure similar to that of 
R2 both of them have a transient interval also. We refer to these two transient bands as 
peripheral transients since in addition there is a central transient in the interval X^^c- In 
general, as discussed below (@y), at the n th bsp R n has 2 n_1 — 1 transient bands, of which 
2^—1 _ 2 are peripheral transient bands. At each bsp all transient bands except the central 
one are rescaled versions of the central transient bands at previous band-splitting points. 

Under the map R n the inverse image of any point in the central interval X n Q is 
contained within the interval itself. This implies that if a function initially has no support 
in X nt c, it will continue to have no support in X n ^Q under repeated iterations by U-& n . 
The spectrum and the form of the eigen/Jordan vector in the central interval at all the 
band-splitting points can be obtained by inspection. We notice that 



where Ht 



(T n (T n _l...(Tl 




(75) 



UR n - X* n ) 3 Xn,C 




} , (76) 
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n=3 



n=2 



n=l 



Figure 7: Band structure at the 1 st , 2 nd and 3 rd band-splitting points. The central 
transient -X^c at the 2 nd bsp transforms into two peripheral transients at the 3 rd 
bsp. 

where Xn,c is the indicator function on X n> c and Xn,b = 1 if ac G [Kn\x c ),^n\x c )) and 
otherwise. The form of the eigen/ Jordan vector in X nj c associated with decay out of 
the central interval is thus (x — The complete form will be determined below. The 

eigenvalues associated with decay out of the central transient at the n th bsp are seen by 



Next we obtain the complete spectrum at the n bsp associated with the decay out of 
all the transient regions onto Q n . This is done by transforming all the central transients at 
band-splitting points of order less than n. Of all the peripheral transient bands at the n th 
bsp, 2 n ~ 2 of these are (n — 2) times rescaled versions of the central transient at R2, 2 n ~ 3 of 
these are (n — 3) times rescaled versions of the central transient at R3 and so on up to 2 of 
these peripheral transients being rescaled versions of the central transient at the (n — l) th 
bsp. This is shown in Figure 7 up to n = 3. We denote the transient eigenvalues at the 
n th bsp by 4>l'\, where I indicates that it was obtained from the central transient at the 



© to be 




(77) 



2G 



(n — Z) th bsp, where I = 0, 1, 2, . . . , n — 2, fc is an integer from 1 to 2', j = 0, 1, 2, 3, ... , and 



as before the right state \<t>^\) is piecewise-polynomial of degree j in each of the pieces. 
The eigenvalues are obtained in a similar fashion to (pH) as 



(78) 



where the 0"'° are given by ([77]). For even values of j there are degeneracies in the 
spectrum while for odd values of j there are no degeneracies. We consider first the states 
associated with the degenerate eigenvalues. 

3.2.1 Transient right states of even degree 



As seen in (78), for a given I (with n and 2j fixed) there are 2 l distinct eigenvalues indexed 
by k. At each integer step (I) all the eigenvalues of the previous step (I — 1) are present and 
2 l ~ 1 new eigenvalues appear. But identical eigenvalues from different steps have disparate 
k values. It is convenient to rearrange the k index so that degenerate eigenvalues share 
the same Table 2 contains </>^'q for all possible values of I and k, arranged to illustrate 
the reordering in k. 





k = 1 


k = 2 


k=3 k = 4 


k = 5 


k = 6 


k = 7 


k = 8 


h,Q 


f 
















<P 


-<P 












<j>k,0 


<P 


-<P 


iip —iip 










j.5,3 


f 


-ip 


iip —iip 


</?exp(f ) 


</>exp(^P) 


(^exp(^) 


(^exp(^) 



Table 2: Transient spectrum at the 5 th bsp with j = 0. The constant ip = 1/05. 



Next we ask if there are independent eigenvectors associated with the degenerate eigen- 
values. Using the same procedure used to find the right eigenvectors decaying out of region 

1 This can be accomplished by choosing the new fc's asssociated with the I th step as k l 2i _ 1 = k 1 ^ 1 and 
k\i = k\i_x + 2 l ~ 1 , where i — 1, . . . , 2 !_1 denotes the order of its appearance and fc? = 1. 



27 



I it can be shown that in the rescaled map, for n > 2 there is no eigenpolynomial of even 
degree with support in the central transient. This means that the right states associated 
with 4>i'2j are Jordan vectors, for n > 2. For n = 2, decay out of the central transient 
is described by eigenvectors for all values of j. Since the right states associated with 
4%j, 4*k' 2 j'- • • ^k'J -2 are transformed versions of these are also eigenvectors. For 

n > 3 decay out of the central transient is described by Jordan vectors for even values of 
j. Therefore all the peripheral transients which are related by conjugacies (|44] ) and (|46| ) 
to these central transients are also Jordan vectors. Hence associated with the transient 
spectrum are Jordan blocks whose sizes correspond to the algebraic multiplicity of the 
eigenvalues. Thus from Table 2 we see that at the 5 th bsp there is a 4 x 4 Jordan block 
'2j, a 3 x 3 Jordan block associated with $2 2j 

i „-,,i, jJ , 2 . 

n—lxn—1. Since 0^' o is the largest eigenvalue with modulus less than 1 in the spectrum 
of Un n , it corresponds to the slowest decay mode. Since there is a Jordan block associated 
with </>i'q this decay is modified exponential (polynomial factors in t times exponential 
decay) . 

In general (for even j) we have Jordan vectors for I < n — 3 and eigenvectors for 
I = n — 2 as 

^Rnl^fei') = d'lCi?) + ICS 1 ) n>3, /<n-3 

OrJ^S"*) = 4>Zf\^lf) n ^ 2 - ( 79 ) 

By inspection we have already obtained the form of the Jordan state in the central region 



5 5 1 tVi 

associated with 4>\2ji a 3 x 3 Jordan block associated with </>22j an< ^ so on - ™ ^ ne n ^ S P 
the largest Jordan block in the transient spectrum is associated with ^"'^ and is of size 



in (|76|). We have 

\4>x%) = a„, 2 j (X - X* n f J Xn,C + fn,2j(x) (80) 

where f n ,2j{x) is piecewise polynomial with degree 2j over each of the S n intervals (ex- 
cluding the central interval X n ^) at the n th bsp. Since a Jordan vector multiplied by a 
scalar doesn't remain a Jordan vector (with respect to the same eigenvector in the block) 
we do not have the freedom to choose the a n ,2j's to be 1 as we did in (^6|). Using (|7(f ) and 



2S 



( |79| ) we know that if (80) is to be a Jordan vector the arbitrary functions f n ,2j(x) must 
satisfy 



,71,1 



UR n f n ,2j(x) = 4>i§j [fn,2j(x) ~ a n fr (x - X* n f 3 Xn,b\ + \4>1%) ( 81 ) 

A formal approach to determine the a n .2/s and f n ,2j(x) is described in Appendix D. 

To find the eigenstates corresponding to the peripheral transients we transform the 
central transient eigenstates at the 2 nd bsp to all the higher band-splitting points. We 
define 

b^ 2 = (C 2 T 2 ) 1§ - < 82 ) 

and n o - n( T n _ 1 ...o-3 is an ordered product of n — 2 operators Aj and o^^Ji (i = 3,4, . . . , n). 
If Oi = the i th place in the product is taken by Aj and if Gi = 1 the i th place in the 
product is taken by B^'^ - /. Following the same procedure used to obtain (|6l|) we find 
that the transient eigenvectors, \<P2'2j~ 2 )^ ( n — are given by 

IC 2 T 2 ) = En^-,...3l^°) (83) 

M 

To find the Jordan states corresponding to the peripheral transients, \<Pk'2j) f° r 1 — ^ — 
n — 3, we use an approach similar to that used to find the Jordan vectors describing decay 
out of the central transient, l^"'^)- To clarify the notation we note that |^' 2 j) does not 
have support in the central transient band, but has support over all the other transient and 
attracting bands. Similarly |0^2j) has support on all bands except the central transient 
band and the two transient bands that are related by a single transformation to the central 
transient at the (n — l) th bsp. This pattern continues and |</ , 5?' 2 j~ 2 ) has support on all the 
attracting bands and the 2 n ~ 2 transient bands that are related by transformations to the 
central transient at the 2 nd bsp. Following (|80| ) and using the fact that the peripheral 
transients at the n th bsp are transformed versions of the transients at the (n — l) th bsp 
we write 

Wkij) = a k,2j (^n\^k/2l2j) + $kifin\<t>"kW~2j)) + <2j( 2; ) ( 84 ) 

where g^ ,l 2 j(x) is piecewise polynomial of degree 2j. The polynomial g^ ,l 2 j(x) can be ob- 
tained using the procedure similar to the one described in Appendix D. We present here 



29 



the solution only for the simplest case 



n,n—3 
'k,2j 



( 



) 



2 6 



n,n-3 Wk'Pj 
'k,2j 



\,n,n-2 



(85) 



where k' = k + 1 if k is odd and k' = k — 1 if k is even. The recursion relations for the 
other Jordan vectors involve more terms, and in general to find all the Jordan vectors at 
the n th bsp one must know all the Jordan vectors at the (n — l) th bsp. Also at a particular 



can be obtained directly from (g^) without regard to this order. 
3.2.2 Transient Right states of odd degree 

For odd values of j the transient spectrum at the n th bsp, given by equation (^), is 
nondegenerate. We first find the eigenvectors with support in the central transient bands, 
at all the band-splitting points. By inspection we have already obtained the form of the 
eigenvector in the central transient in ( |76| ) so that 



where f n ^j+i{x) is a polynomial of degree 2j + 1 in each of the S n intervals excluding the 
central interval. For this expression is to be an eigenvector f n ^j+i{x) must satisfy 




101,'jy+l) =(x~ <) 3+ Xn,c + fn,2j+l{x) 



(86) 



UR n f n ,2j+l(x) = /„,2j+l(x) - (X - X* n ) 2j+1 Xn,b ■ (87) 



This equation for f nj 2j+i{x) is similar to equation (|8l]), except that it has one term less 
on the rhs, and may be solved using the procedure outlined in Appendix D. 



Once the l^i^j+i) are known the eigenvector corresponding to any |<^' 2 j + i) can be 
written as 






(89) 
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If <7j = the i place in the product is taken by Aj and if <7, = 1 the i place in the 
product is taken by ,l 2 - f . The sum is over all possible strings of O's and l's of length I. 

3.2.3 Left States 

The left states (4% '.| form a bi-orthonormal set with all the previously obtained right states 
at the n th bsp as 

(€jl A fcV> = ( 90a ) 
(^kjl^k'f) = fill'fikk'Sjj'- (90b) 

Among the right states only the states with I = have support in the central transient 
band. From (|^) we see that the associated left states are 

<#&l = —^>(s (91a) 

(<t>i$ j+ i\ = -* {2J+1) (* ~ <) • (91b) 

The left states corresponding to the peripheral transients at the n th bsp, are found 

by transforming the left states at the (n — l) th bsp, (^n^l ~j^\ m Iben the entire set at the 
n th bsp has to be orthonormalized using a Gram-Schmidt procedure. 

3.3 Back to T n (x) 

Going back to the tent map T n (x), we transform all the right states of R n (^) by U,-i and 
the left states by K^. The states describing decay out of [T n {x c ), 1] are null states. The 
eigenvalues describing decay out of [0, (x c )] are = f^V and have associated 
polynomial eigenvectors of degree j. For j even this part of the spectrum overlaps with 
the spectrum describing decay of transients of R n . There are Jordan vectors associated 
with this part of the spectrum; denoting them as \(j)2j) we have 

DtJ^>=^|^> + |^&> (92a) 
UT n \^ j+1 ) = ^ j+1 \^ j+1 ). (92b) 
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These right states can be determined by an extension of the methods used in Section 2.2 
to determine the eigenstates describing decay out of interval I for the tent map at the first 
bsp. 

4 Conclusion 

We have presented the generalized spectral decomoposition of the Frobenius-Perron op- 
erator of the tent map at the band-splitting points. The right eigenstates are piecewise- 
polynomial functions and the left eigenstates are generalized functions. The spectrum is 
discrete and gives the characteristic decay times of the map. From the decomposition 
one can calculate correlations of arbitrary polynomials (as well as functions expandable 
in terms of the polynomial eigenstates). Furthermore, since the modes corresponding to 
transient decay onto the attractor have been obtained, the full nonequilibrium dynamics 
of initial probability densities is accessible. 

The slowest decay mode, corresponding to the eigenvalue a" 1 at the n th bsp, describes 
decay onto the attractor. At the n th bsp there is an n x n Jordan block associated with this 
eigenvalue and therefore the decay is modified exponential. The asymptotic periodicity of 
the map is clearly reflected in the spectrum as at the n th bsp, all the n th roots of unity 
are part of the spectrum. Our analytic solution of density evolution in this system may 
be useful for comparision with the behavior of systems governed by the Ginzburg-Landau 
equation since a component of its dynamics Q can be reduced to the tent map. 
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Appendix A: Topological conjugacy 

In this appendix we review the spectral decompositions of maps related by a coordinate 
transformation Q. Let T : X — > X be a map defined on the interval X. Transforming 
the interval X by the one-to-one, onto, continuous function : X — ► Y gives a new map, 
S : Y — » y. This map is determined as 

Vt+i = 4>(xt+i) = 0(T(x t )) = S(yt). (93) 
Using that has an inverse gives 

0(T(x*)) = 4>{T{4>- \4>{x t )))) = 4> o T o ^r 1 ^), (94) 

so that 

S = ^oTof 1 , (95) 

The maps T and S are said to be topologically conjugate to each other. 
The Koopman operator, Kg, corresponding to S is given from (|95|) by 

K s = K^KtK^ (96) 

where we have used the fact that K^-i = K^ 1 . The Frobenius-Perron operator, Us, 
coresponding to S is the adjoint of K$. Taking the adjoint of ( p6| ) and using (K^ 1 )^ = 
{K\)~ l gives 

U s = U^U T U^. (97) 

Since Us and Ut are related by the similarity ( p7| ) the spectrum of Us is identical to that 
of Ut and eigenstates transform as 

|A n ) s = U<f, |A n ) T , (98) 

where we use a Dirac-style bra-ket notation for the states. From (|9^) the left states 
transform as 

(A n | s = K^-i (X n \ T . (99) 

Jordan states of the maps are also related as in (|9^) and (^) and both the algebraic and 
geometric multiplicities of the eigenvalues are preserved under conjugacy. 
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Appendix B: The tent map with unit height 



The Frobenius-Perron operator of the tent map with unit height is given by 



Ut p(x) 



+ P 



x 



(100) 



The operator Ut admits polynomial eigenstates with support on the whole unit interval. 
Associated with eigenpolynomials of order 2j are the nonzero eigenvalues 2 _2j . There is 
an infinite degeneracy of the eigenvalue with an independent odd-order eigenpolynomial 
associated with each occurence of the eigenvalue. Thus, the odd-order eigenpolynomi- 
als are not unique but we choose them as Euler polynomials so that the associated left 
eigendistributions take a simple form. The right eigenvectors of Ut are 0> |l2| 

|2-^')t = B 2j (x/2) (101a) 
|0 2i+ i) To = E 2j+1 (x), (101b) 

where Bj(x) is the Bernoulli polynomial of order j and Ej(x) is the Euler polynomial of 
order j [O]. The corresponding left states are 



(2- 2 ^| To = 2^B 2j {x) 



(0 2 j+i| To 
where Bq(x) = 1 and for j > 1 
B 2j (x) = 



1 



(2j + l)! 



5^ +1 \x-l), 



1)2,-1 



x-1) 



c(2?-l) 



(X) 



(102a) 
(102b) 

(103) 



where the action of 5 < ^ l \x — c) on a sufficently differentiable function f(x) is given by 



/ dx5 { r ) (x-c)f(x) = ]im(-irf 



m Am) 



c±e), 



(104) 



for a < c < b and e is a positive infinitesimal. 

The time evolution of a density is expressed in terms of the spectral decomposition of 
Ut as 



£4 fa) = E [(^ 

3=0 



)(2- 2 i\p) + S t>0 \0 2j+1 )(0 2j+1 \p) 



(105) 
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where the bilinear form is defined by 



</|<7>= f 1 dxf*(x)g(x). (106) 
Jo 

Appendix C: Calculation of transient right states 

To determine the functions gnj(x) and ffm : j(x), which appear on the rhs of (|37|) we expand 
9n,j(x)xn + 9ni,j{x)xui m terms of the eigenstates given in Table 1 of Ut 1 on the attractor 
as 

Li/2J b'/2j 
9n,j(x)xn + 9m,j(x)xm = E a i,j\+ 2 ~ l )n + E b ij\~ 2 ~ l )n 

i=l i=l 



+ E c i j|0ai+ 1 >n+ E di tj \0j 2i+1 ) n , (107) 



i=i i=i 

where [^J denotes the integer part (floor) of the real number x. Then acting with Ut 1 
gives 

Li/2J b'/2j , 

i=l i=l 

+ E *j|02i+i>n. (108) 

i=l 

We substitute fll07| ) and ( |lQ8|) into (§?]) and act on @ by all the left states on the 
attractor in succession. Using orthonormality we obtain the following equations for the 
expansion coefficients: 



a i,j 



+ a iyj = 2- (j+1 ^ 2 aij (109a) 



b f ~ hi = -2- (j+1)/ X- (109b) 
dij + Tij = 2-U +i y 2 c itj (109c) 
djj = 0, (109d) 
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where 

anj = 2-^' +1 )/ 2 (+2- i |x J xn) (HOa) 

Aj = 2-^' +1 )/ 2 <-2- < |^xn>=Oi J - (HOb) 

7ij = 2-^'+ 1 )/ 2 (0 i |^ X ii) ; (HOc) 

and we used (Ojja^xn) = because (OjJ has support only in xiii- Explicit evaluation of 
(+^\x 3 Xn) gives 



Xii) = < 



Evaluation of (0j|a; J xil) gives 



5R 4 5 - ;y ((2-^-l) i = 

2t - 1 > j (111) 

(2i)!(j-2i+l)! ^ Z i; Z« 1 <. J. 



<0 i \x'xn) = { ^y-^- 1 ^ 1 ~ J (112) 

otherwise. 

These results are then used in (109) to determine the expansion coefficients for the transient 
eigenstates with support in region I. 

Appendix D: Transient right states at higher bsps 

We expand the arbitrary functions f n ,2j(x) in terms of the transient and non-transient 
eigenvectors of Uji n . The non-transient eigenvectors are given in ( |6l| ) and ([7^) . The tran- 
sient eigenvectors will be transformed versions of the central eigenvectors at the previous 
band-splitting points. The expansion is 

j 2 n j 1 

fnxW = EE ®^> + EEi?iW 
k=ij' =1 k=if = i 

n-2 2 l j n-2 2 l j-l 

+ E E E <?Kj) + E E E <f K^- ^ 
i=i k=ij'=i i=i fc=ij'=i 
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Since d^'j 2 ' 2 -' is the coefficent of the eigenvector of the Jordan block it can be set to zero. 
Applying Un n to the function f n ,2j(x) we get 

2" j 2 n — 1 j—1 

unjnM*) = E E ^, 2 /l A *,2/> + E E <5S/I<W 

fc=l/=l k=l j'=i 

n-2 2' j n-2 2 i ~ 1 j 

+ EEE ^TOO +EEE 

i=l fc=lj' = i Z=2 fc=l j' = i 

n-2 2' J 

i=i k=if = i 

We substitute ( |114| ) and (113) into ( |3lD and hit both sides of the equation with (A" .,\, 



<°fc,jy'+il' ^fci'+i 1 and ^fei' 1 successivel y- Lettin S 



"J? = ( A fc,2j'l (^-^nfiXb.n) (H5a) 
= (0k,2? + l\(x-O 2j Xb,n) (H5b) 
7^ = (<P n k i f \ - x* n f Xb ,n) (H5c) 
we obtain the following equations for the expansion coefficents a n j, b™' 2 j , c™' 2 /, cC' l fi and 

' fcjj 

n,i,2j 



„n,2j _ ,n,0 (n,2j _ R n,2j 

c k+i,j' ~ rij {.%? p k,f 



in,0 n,l,2j i 
a n,2j<t> h2 jJl,j = 1 



n,l,2j,n,l 
c k,j' yk,2j'+l 



The equations for d?''', 2 "' differ depending on the values of /,£; and f. For k = 1, f = j 



and I = 2, 3, . . . , n — 2 



When j 7^ j', A; 7^ 1 and / = 1 we have 



jn, I— 1.27 j.n.0 n,l,2j 

<hj = a n,2j<Px, 2 j1lj 

d n,n-2,2j = Q (m) 



jn.l,2j / /n,0 iTi.l \ ^">0 n,l,2j /-. -. d \ 
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When j ^ f , k ^ 1 and I = 2, 3, . . . , n — 2 we have 

%j> \<Pl,2j ~ <Pk,2j> ) ~ d k,j> - a n,23<Pl,2j%j> ■ l iiy i 



The equations ( |116| ) - (119) for a n ^j, b^'p, c£'p d 1 ^ 1 -^ and e^'^', 2 "' are either uncoupled 



or are coupled in a simple manner and can be solved explicitly to find the expansion 
coeffients. Plugging these coefficents into ( |1 13|) and (^) we get all the Jordan states with 



support in the central transient, for all n > 2. For n = 2 set all the d l ^ k ., = 0, e^.', 2jf = 



and d2j = 1 and solve for \? n k •,, fcjV from ( 119 ) to obtain the eigenvectors with support 
in the central transient at the second bsp. 

References 

[1] A. Lasota and M.C. Mackey, Chaos, Fractals and Noise (Springer- Verlag, New York, 
1994). 

[2] H.H. Hasegawa and W.C. Saphir, Phys. Rev. A, 46, 7401 (1992); P. Gaspard, J. 
Phys. A, 25, L483 (1992); I. Antoniou and S. Tasaki, Physica A, 190, 303 (1992); 
I. Antoniou and S. Tasaki, J. Phys. A, 26,73 (1993). 

[3] H.H. Hasegawa and D.J. Driebe, Phys. Rev. E 50, 1781 (1994); P. Gaspard, Chaos, 

3, 427 (1993); S. Tasaki, I. Antoniou, and Z. Suchanecki, Chaos, Solitons & Fractals, 

4, 227 (1993). 

[4] D.J. Driebe Fully Chaotic Maps and Broken Time Symmetry (Kluwer Academic Pub- 
lishers, Dordrecht, to appear). 

[5] H.G. Schuster Deterministic Chaos (VCH Publishers, New York, 1988). 

[6] H.T. Moon, Phys. Rev. E, 47, 772 (1993); H.T. Moon and M.J. Kim, Int. J. Bifurc. 
Chaos, 7, 1141 (1997). 



38 



[7] S.J. Shenker and L.P. Kadanoff, J. Phys. A, 14, L23 (1981); T. Yoshida, Prog. Theo. 
Phys., 72, 895 (1984); R. Brown, C. Grebogi and E. Ott, Phys. Rev. A, 34, 2248 
(1986). 

[8] H. Shigematsu, H. Mori, T. Yoshida and H. Okamoto, J. Stat. Phys., 30, 649 (1983); 
T. Yoshida, H. Mori and H. Shigematsu, J. Stat. Phys., 31, 279 (1983). 

[9] M. Dorfle, J. Stat. Phys., 40, 93 (1985). 

[10] N. Provotas and M.C. Mackey, Physica D, 53, 295 (1991). 

[11] D.J. Driebe and G.E. Ordonez, Phys. Lett. A, 211, 204 (1996); G.E. Ordonez and 
D.J. Driebe, J. Stat. Phys., 84, 269 (1996). 

[12] R.F. Fox, Chaos, 5, 619 (1995); I. Antoniou and B. Qiao, Phys. Lett. A, 215, 280 
(1996). 

[13] J. Heidel, Phys. Lett. A, 143, 195 (1990). 

[14] M. Abramowitz and LA. Stegun, eds., Handbook of Mathematical Functions, (Dover 
Publications, New York, 1972). 



39 



